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Abstract 

Nonlinear hydrodynamics is used to evaluate disorder-induced corrections to the 
vortex hquid tilt modulus for finite screening length and arbitrary disorder geometry. 
Explicit results for aligned columnar defects yield a criterion for locating the Bose 
glass transition line at all fields. 
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The vortex phase diagram of high temperature superconductors shows a rich 
diversity of phases, including vortex lattices, liquids and glasses. Novel types 
of glasses are also possible because of pinning in disordered samples (for a re- 
view, see Refs. [1-3]). Much progress in the theoretical understanding of vortex 
behavior at low fields has been made by employing the formal analogy of the 
statistical mechanics of (2 + l)-dimensional directed lines with the quantum 
mechanics of 2d bosons. [4] In this mapping, the vortex lines traversing the 
sample along the direction z of the external field, Hq = zHo, correspond to 
the imaginary-time world lines of the quantum particles. [4] The thickness, L, 
of the superconducting sample is the inverse temperature, Pb^, of the bosons, 
the vortex line tension, ei, represents the boson mass, m, and thermal fluctu- 
ations in the vortex state oc fc^T map onto quantum fluctuations oc ^ in the 
boson system. This "boson mapping" has been particularly useful for under- 
standing the properties of vortices in superconducting samples with aligned 
damage tracks from heavy-ion irradiation. This type of disorder yields a low- 
temperature Bose-glass phase where every vortex is trapped on a columnar 
defect [5,6] and vortex pinning is strongly enhanced. [7] The transition at Tbg 
from the entangled vortex liquid to the Bose-glass is continuous and is signalled 
by the vanishing of the linear resistivity and of the inverse tilt modulus, C44 . 
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The drawback of the boson mapping in the form used by Nelson and cowork- 
ers [4,6,8] is that intervortex interactions are assumed to be strictly local in 
z, the magnetic field direction. This restricts the application of the results to 
low fields. It also renders problematic the evaluation of the tilt modulus C44, 
which measures the linear response to a transverse tilting field normal to 
Hq. The long-wavclcngth tilt modulus of the vortex liquid can be written as 
C44 = C44 + B^/(A7i), where C44 = noii, with uq = B / (j)Q the vortex density, 
is the single vortex part and the second term represents a compressive con- 
tribution. By assuming local interactions along the compressive part of C44, 
which dominates at high fields, is neglected entirely. With this approximation, 
the mathematical analogy between vortices and 2d bosons can be exploited 
further to show that the inverse tilt modulus maps onto the superfiuid density, 
of the bosons, with 
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The transition as T ^ from the entangled vortex liquid to the Bose 

glass corresponds then to the transition from a boson superfiuid to a localized 
normal phase of bosons. Tauber and Nelson have evaluated perturbatively 
the reduction of Ug due to various types of disorder. [8] They also showed that 
for aligned columnar disorder such a perturbative calculation yields a useful 
criterion for locating the transition line Tbg{B) at low fields. 

Larkin and Vinokur [9] argued that a generalization of Eq. (1) that incorpo- 
rates the compressive part of C44 can be obtained by using the nonlocal map- 
ping on vortex lines onto 2d bosons introduced some time ago by Feigel'man 
and collaborators. [10] These authors showed that the fully nonlocal London 
model of interacting vortex lines can be mapped to a system of 2d charged 
bosons coupled to a massive photon field. The duahty between vortices and 
bosons translates into 



where Ug is defined here in terms of the polarization function of the fictitious 
gauge field. Near a Bose glass transition where — > the second term of Eq. 
(2) dominates and results in the divergence of the tilt modulus. Feigel'man et 
al. carried out a perturbative calculation of the reduction of from intervortex 
interactions in a clean material (in the limit A 00) [10], but the calculation 
of the reduction of Ug from disorder in the context of the charged boson model 
is cumbersome and does not provide much physical insight. 

An alternative approach for modeling interacting vortex arrays is hydrody- 
namics, which has proved very useful to describe the long wavelength proper- 
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ties of flux-line liquids at high fields. [11] Hydrodynamics provides a physically 
transparent formulation that naturally incorporates the nonlocality of the in- 
tervortex interaction. The goal of this note is to show how hydrodynamics can 
be used in a transparent way to evaluate disorder-induced corrections to the 
wave- vector dependent tilt modulus for finite values of the screening length A 
and arbitrary disorder geometry. The connection of the hydrodynamic formu- 
lation to the charged boson formalism will also be discussed. Explicit results 
are presented for aligned columnar defects. 

The hydrodynamic free energy of the flux-line liquid is a functional of two 
coarse-grained flelds, the local areal density of vortices, n{r), and the tilt field, 
t(r), which measures the local deviation of a volume of fiux liquid from the 
direction of the external field, Ho.[12] It is given hj F = Fq + Fd, where Fq 
is the free energy in the absence of disorder and Fd describes the coupling to 
quenched defects. The disorder-free contribution is 

= r-Wf^i^l K4(q)|t(q)r + c?,(q)|Mq)r} , (3) 

r ^ ' ^ q 

where 5n(r) = n(r) — no. The density and tilt field are related by the familiar 
constraint that flux lines cannot start nor stop inside the sample, 

a,n(r) + V±-t(r) = 0. (4) 

The bare elastic constants, c'f^ (the compressive part of the tilt modulus) 
and (the bare compressional modulus), are determined by the intervor- 
tex interaction. For isotropic materials they are simply cf^iq) = c^'J(q) ~ 
(S^/47r)/(l -I- q^X^), with A the screening length. [13] Quenched disorder from 
material defects couples to the flux-line density and gives a contribution 

Fd = J VD{r)Sn{r). (5) 

r 

The random potential Voir) is taken to be Gaussian, statistically homoge- 
neous, and isotropic in the xy plane so that 

SVDici)5VDici') = A(gx, g.)r2()q,+q'„o , (6) 



where SVD{r) = VD(r) — Voir), and the overbar represents the disorder aver- 
age. The correlator A.{q±, q^) depends on the geometry of disorder and will be 
specifled below for the case of interest. 

The hydrodynamic free energy Fq goes beyond the Gaussian hydrodynamic 
model commonly used in the vortex literature [11] as the flrst term on the right 
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hand side of Eq. (3), describing the "kinetic energy" part of the vortex inter- 
action, incorporates non-Gaussian terms. In a recent pubhcation we showed 
that this non-Gaussian hydrodynamics is precisely equivalent to the charged 
boson model of Feigel'man and coworkers. [10] Such a "nonlinear hydrodynam- 
ics" was used in Ref. [12] to evaluate perturbatively the enhancement of C44 
from interactions in a clean material. It was shown there that nonlocality is 
crucial to yield a correction to C44 that remains finite for L 00. The same 
model is used here to evaluate corrections to C44 from disorder. In particular, 
for aligned columnar defects our calculation yields a criterion for locating the 
transition line Tbg{B) at all fields. 

The tilt modulus, C44, can be expressed in terms of the tilt field autocorrelation 
function, Tjj(q) = (tj(q)tj(— q)), where the brackets denote a thermal average 
with weight ~ exp[— (Fq-I-Fd) //cbT] , to be carried out subject to the constraint 
(4). It is given by 



with P^^(qx) = 5ij — q±iq±j the familiar transverse projection operator, and 
q± — (l±/q±. When only quadratic terms in the fluctuations 5n and t are 
retained, the resulting Gaussian free energy, Fq, is Fq = Fqg + F^,, where 
Fqg is obtained from Eq. (3) by the replacement |tp/n(r) |tp/no. The first 
two terms on the RHS of Eq. (3) can then be combined to define a bare tilt 
modulus C44 = nil + C44 and Eq. (7) is simply an identity. To Gaussian order 
there is no coupling between the density field and the transverse part of the tilt 
field that determines C44. As a result, disorder, which couples to the density, 
docs not change the tilt modulus. Non- Gaussian terms in the hydrodynamic 
free energy of Eq. (3) do introduce a coupling between transverse tilt and 
density, yielding a renormalization of the tilt modulus. Specifically, the free 
energy is written as F — Fg + 5F, with 



nlksT 



hm hm P^(q^)r,,(q) , 



(7) 





By treating the non-Gaussian part 6F perturbatively, we obtain 



C44 0^4 L clJnlnQy 



(9) 



with 
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The brackets {■■■)g in Eq. (10) denote a thermal average with weight ~ 
cxp[— Fg/A;bT], subject to the constraint (4), followed by the average over 
quenched disorder. The correction has been denoted by n„ because it has the 
suggestive interpretation of a normal-fluid density. Notice that its expression 
is formally identical to that obtained in [12] in the absence of disorder. Disor- 
der only enters here in the expressions for the Gaussian correlators, which can 
be found for instance in Ref. [12]. A more general expression for the renor- 
malized tilt modulus at finite wave vector can also be obtained by the same 
methods. [14] In the dilute limit, X « a, with a = 1/y^ the intervortex spac- 
ing, the 2;-nonlocality becomes insignificant and Eq. (10) reduces to the result 
obtained for instance by Tauber and Nelson using the local boson mapping. [8] 

We now focus on the case of correlated disorder created by heavy ion ir- 
radiation yielding rectilinear damage tracks aligned with the field Hq. This 
corresponds to A(q) = AL6q^fi, where A Ugb'^/cP, with Uq the depth of 
the pinning potential well, b its radius, and d the average distance between 
columnar defects. [6] In this case the disorder contribution to n„ from disorder 
takes a particularly simple form, given by 



where G{q±) = 1 -|- c^4^{q±, Qz = 0)/cl4^{q±, — 0), and eB{q±) corresponds to 
the spectrum of 20? bosons with screened interactions. 



with V{q±) — Vb/(1 -|- X']_q'^) the screened boson interaction and Vq — (t>i/^T^- 

The function G{qi_) depends only weakly on g_L, varying between 1 and 2, 
and is the only manifestation of the 2;-nonlocality of the intervortex interac- 
tion. If we let G{q±) = 1 and negelct the screening, i.e., V{q±) = Vq, our 
result becomes identical to the disorder- induced renormalization of I/C44 ob- 
tained by others. [15,8] As expected, for the case of aligned columnar defects 
the 2-nonlocality is not very important. The screening of the boson interac- 
tion incorporated in V{q±) is, however, important at high vortex densities. 
This becomes clear by plotting the boson spectrum as a function of wave 
vector, as shown in Fig. 1. At high density the screening of the interaction 
yields a flat region in the spectrum. The wave vector integral in Eq. (11) is 
dominated hj q± ~ ksz- At low density the main contribution to the in- 
tegral comes therefore from a region where the spectrum is phonon-like, i.e., 

^b{(1±) I {ksT) ~ (l±\/nQV{q±, T)/ei, as obtained from the theory of uncharged 
boson superfluids.[16] At high density, however, the main contribution to the 
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integral comes from a region of wave vectors where the spectrum is plasmon- 
hke, i.e., eB{q±) / {kBT) ~ i/iTmo, as appropriate for a charged supcrfluid. This 
was actuaUy recognized earher by Larkin and Vinokur, who in order to use the 
results of the local boson theory at high fields proposed an ad hoc formula for 
the Bogoliubov spectrum that captures the dense liquid physics by interpolat- 
ing between these two limits. [9] Hydrodynamics naturally provides a simple 
and unified description of the behavior of fiux-linc liquids in the presence of 
disorder that applies over a wide range of densities. 

As discussed in Ref. [12], the perturbation theory breaks down at high fields 
and temperatures [kBT/iia > 1). An approximate expression for the tilt 
moduls of the form proposed by Larkin and Vinokur [9] can be obtained 
by treating the nonlinearities in a mean field approximation, as done by 
Feigel'man and collaborators for the charged boson superfiuid.[10] The tilt- 
tilt autocorrelation function is evaluated by retaining only transverse fields 
and using a Hartree-type approximation for the non-Gaussian terms. [14] This 
yields 



C44 cfj no + 61/ (no - Un) ' 



with Un given approximately by Eq. (10). 

The Bose glass transition line Bbg{T) can now be evaluated as the locus of 
points in the {B,T) plane where n^/riQ = 1, corresponding to C44 — >• 00. In 
general the integral in Eq. (11) has to be evaluated numerically. The result- 
ing Bbg{T) phase line is shown in Fig. 2. Analytical results can be obtained 
in the limit of low and high field. For a dilute liquid (A >> a), the screen- 
ing of the interaction can be neglected and the Bogoliubov spectrum can be 
approximated as linear in q±. The integral can then be evaluated and gives 

, , ln(l/nA2) /r*\4 /Tc-T\^ 
™«B,^^(-) ~(^), (14) 



where = (j)Q/cP is the matching field and ksT* = b^eiUo)^^"^ is a char- 
acteristic pinning energy scale, in agreement with earlier work by other au- 
thors. [15,9,8] In dense liquids (A >> a), the Boguliobov spectrum is approxi- 
mately gf_L-independent. Again, the transition line can be obtained analytically, 
with 



which agrees with the earlier result of Larkin and Vinokur [9]. 
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Hydrodynamics provides a physically transparent framework for evaluating 
disorder-induced corrections to the elastic constants of the vortex liquid. It 
naturally incorporates the full nonlocality of the intervortex interaction and it 
allows us to compute finite- wavevector elastic constants for arbitrary disorder 
geometry. In particular, the nonlocality in the field (z) direction is expected 
to be important for splayed columnar defects, where terms coupling disorder 
directly to the tilt field may need to be incorporated in the coarse-grained 
theory. This work will be presented elsewhere. [14] 
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Figure Captions 

Fig. 1. The boson spectrum e* = [eB{Q±) /{kBTkBz)]'^ versus wavevector 
K = q_L/kBz for noA^ = 0.1 (solid line) and tt-qA^ = 10. (dashed line). 



Fig. 2. The Bose glass transition line as obtained from the condition n^/no = 1. 
The vertical axis is B* = Bbg{T)/B^ and the horizontal axis is t = T/Tc. The 
parameter values used correspond to B^ = 2.3T, b = 20A, and T* = O.lTc at 
T = 0. The transition line is well approximated by Bbg ~ as in Eq. (15) for 

all but extremely low fields. The dashed line represents the the mean field Hc2. 
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